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Abstract. We study the chiral phase transition by the non-equilibrium propagation of the 
sigma field. A quark fluid acts as a heat bath in local thermal equilibrium and evolves 
fluid dynamically. We allow for dissipative processes and fluctuations since the sigma field 
is propagated according to a Langevin-type equation of motion. Non-equilibrium fluctuations 
at the first order phase transition lead to an increase in the intensity of sigma excitations. 



1. Introduction 

The QCD phase transition is subject to intensive theoretical and experimental research. Besides 
theoretical indications of a QCD critical point [TJ [2] experimental signatures are proposed to 
discover the suggested critical point in experiment [3j 2] . The fast dynamics during the expansion 
of the system created in a collision makes it important to consider non-equilibrium phenomena. 
At a critical point relaxation times become infinite. A system that cools rapidly through a second 
order phase transition is necessarily driven out of equilibrium. Any signal for a critical point is 
crucially weakened [5] . A non-equilibrium situation, however, leads to interesting phenomena at 
a first order phase transition, such as nucleation [6] and spinodal decomposition [7]. We embed 
a field theoretical model with a phase transition into a realistic expansion of the hot and dense 
matter and extend similar models [HJ 9j by a complete treatment of fluctuations and dissipation 
of the chiral order parameter coupled to the fluid dynamic expansion. While the chiral field is 
damped energy dissipates into the quark fluid. Random kicks by the noise term also contribute 
to the energy and momentum exchange between the field and fluid. We introduce a stochastic 
source term into the fluid dynamic expansion to conserve energy of the whole system. 

2. Chiral fluid dynamics 

2.1. Semiclassical equation of motion for the sigma field including dissipation and noise terms 
We use the linear sigma model with consituent quarks as a low energy effective model of QCD. 
It exhibits a chiral phase transition at different strengths, a crossover, a critical point and a first 
order phase transition |10j 



C = q [i 7 0„ - g (a + i 75 T7r)] q + \{d^a) + \{d^d^) - U(a, vf) , 



(1) 



with the constituent quark field q = (u,d), the coupling g between the quarks and the chiral 
fields. The potential is given by 

U (a, tt) = ^ (a 2 + vr 2 - v 2 ) 2 - h q a - U , (2) 

where (a) = f n = 93 MeV and (tt) = 0. The explicit symmetry breaking term is h q a = f^m 2 
with m n = 138 MeV. Thus v 2 = f 2 -mir 2 /X 2 . A 2 = 20 yields a sigma mass m 2 = 2A 2 / 2 + m 2 sa 
604 MeV. For zero potential energy in the ground state Uo = m^/(4A 2 ) — f 2 m 2 . 

In the rest of this paper we will concentrate on the order parameter of the chiral phase 
transition, the sigma field and set tt = {tt) = for all times. The equation of motion for the 
sigma field including dissipation and noise is given by 



SU 

Here, the chiral condensate to one-loop level is 



d^a+ — +g(qq) a + r ] d t a = C (3) 
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np(p) is the Fermi-Dirac distribution, d q = 12 the degeneracy factor and E = \/p 2 + g 2 a 2 the 
energy of the quarks. 

Our choice of the damping coefficient r\ and the noise term £ is motivated by the results of 
two- loop calculations in a simplified 0(4) chiral model \Y\\ I12j . In a first approximation we 
assume that r\ is independent of temperature and mass changes and that the noise is white. 
Inspired by [13] we use r\ = 2.2/fm. The dissipation-fluctuation theorem then gives 

2T 

<£(*)) =0 and mat')) = — V d(t - t') . (5) 

2.2. The effective potential and the equation of state 

The pressure of the quark fluid is given by the thermodynamic potential. To one-loop level it is 

V eS (a,«,T) = ~logZ = -2d q T J log(l + e~# ) + U (<r,7f) . (6) 

At fj>B = the strength of the phase transition can be tuned by varying the strength of the 
coupling g. For a small coupling the transition is a crossover, for g = 3.63 the potential becomes 
flat as for a continuous transition. And for g = 5.5 the phase transition is discontinuous. 
The equation of state for the fluid dynamic expansion is obtained from 

p(a,7f,T) = -V eS (a,jt, T) + U{a,jr) (7) 
e{a,ir,T) = T — p(a,vr,T). (8) 

2.3. The stochastic source term 

The equations of relativistic fluid dynamics are 

9,(TZ d + TZ d ) = (9) 
By using the explicit form of the equation of motion ([3]) we derive for the source term 

S" = -dpTZd = ~(9^a + 5 -^)d v a = -(-g(qq) a - rjd t a + fiFa (10) 
The damping and the noise term enter the source term which is therefore also stochastic. 



3. Numerical results 

We solve the fluid dynamic equations ([9]) by a SHASTA code and the field equations © by a 
staggered leap-frog algorithm. The numerical equilibrium solution of equations like ([3j) depends 
on the lattice spacing [T3]. It turns out, however, that the results presented here depend only 
slightly on the lattice spacing, for which we use A = 0.2 fm. The energy density is initiated 
as the equilibrium energy density at T = 160MeV, ellipsoidal in x-y-plane and homogenous 
in z-direction, smoothed by a Woods Saxon type distribution. The sigma field is initiated in 
equilibrium with the quark fluid. All observations shown below are purely due to non-equilibrium 
effects at the phase transition. 

The intensity of the sigma fluctuations is given by 
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(11) 
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The time evolution of is shown in figure ([TJ for a first order phase transition and in figure 
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Figure 1. Time evolution of the intensity 
of sigma fluctuations at a first order phase 
transition. 



200 



150 



average temperature 
average stand dev 
cell temperature 
cell fluctuations 
Tc=123.3MeV 




6 8 10 
time [fm] 



120 



Figure 2. Temperature of the quark fluid 
and standard deviation of the sigma field at a 
first order phase transition. 



© for a critical point. In a first order phase transition the intensity of sigma fluctuations 
between 4 and 11 fm is much larger than at a critical point. This can be better understood 
by looking at the evolution of the temperature and the size of the fluctuations. We show these 
values in a single cell in the middle of the grid and averaged over a volume in the hot region 
of the fireball (V = 0.24 fm 3 ). The standard deviation of the fluctuations of the sigma field in 
V is Ac 2 = 1/V J2v ^ 3 ( cr cq — a ( x )) 2 - F° r a first order phase transition the intensity of sigma 
fluctuations increases around t = 5 fm when the average temperature drops below the critical 
temperature, in figure Here the standard deviation increases to almost 70 MeV. In one 
particular cell in the middle of the hot region one sees that the size of the fluctuation corresponds 
to the difference between the two minima at {ct)t>t c — 10 MeV and (er) V ac = f n = 93 MeV. 
Obviously the fluctuations in a critical point scenario, figure (HJ), are smaller. For a first order 
phase transition we find a second rise in the intensity of sigma fluctuations at around t = 10 fm. 
While the sigma field relaxes from the unstable minimum towards its vacuum value energy is 
transfered to the quark fluid where re-heating occurs. This effect is so strong that the average 
temperature again crosses the critical temperature leading to an increase in the fluctuations. 
Thereafter the sigma field, relaxes slowly towards (o")vac* 
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Figure 3. Time evolution of the intensity of 
sigma fluctuations at a critical point. 
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Figure 4. Temperature of the quark fluid 
and standard deviation of the sigma field at a 
critical point. 



4. Summary 

We presented a dynamic model of phase transitions in heavy ion collisions. We included damping 
and noise originating from the interaction of the chiral fields with the quark fluid. To ensure 
energy and momentum conservation of the coupled system we suggested to introduce a stochastic 
source term into the equations of fluid dynamic. The number of coherently produced soft sigmas 
is much larger at a first order phase transition than compared to a critical point transition. Non- 
equilibrium effects like supercooling and reheating are observed. 

4-1- Acknowledgments 

The authors thank Carsten Greiner, Stephan Leupold, Igor Mishustin and Christoph Herold for 
fruitful discussions and Dirk Rischke for providing the hydro code. This work was supported 
by the Hessian LOEWE initiative Helmholtz International Center for FAIR. M. Nahrgang 
acknowledges finanical support from the Stiftung Polytechnische Gesellschaft Frankfurt. 

References 

[I] Z. Fodor and S. D. Katz, JHEP 0203 (2002) 014 

[2] C. Schmidt [for RBC-Bielefeld Collaboration], J. Phys. G 35 (2008) 104093 

[3] M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81 (1998) 4816 

[4] M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60 (1999) 114028 

[5] B. Berdnikov and K. Rajagopal, Phys. Rev. D 61, 105017 (2000) 

[6] Ph.Chomaz, M.Colonna, J.Randrup Phys. Rev. 389 (2004) 263 

[7] I. N. Mishustin, Phys. Rev. Lett. 82 (1999) 4779 

[8] I. N. Mishustin and O. Scavenius, Phys. Rev. Lett. 83 (1999) 3134 

[9] K. Paech, H. Stoecker and A. Dumitru, Phys. Rev. C 68 (2003) 044907 

[10] O. Scavenius, A. Mocsy, I. N. Mishustin and D. H. Rischke, Phys. Rev. C 64 (2001) 045202 

[II] Z. Xu and C. Greiner, Phys. Rev. D 62 (2000) 036012 
[12] C. Greiner and B. Muller, Phys. Rev. D 55 (1997) 1026 
[13] T. S. Biro and C. Greiner, Phys. Rev. Lett. 79 (1997) 3138 

[14] N. C. Cassol-Seewald, R. L. S. Farias, E. S. Fraga, G. Krein, R. O. Ramos larXiv:071 1.1866 1 [hep-ph] 



